Spin-glass Transition in Bond- disordered Heisenberg Antiferromagnets Coupled with 
Local Lattice Distortions on a Pyrochlore Lattice 
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Motivated by puzzling characteristics of spin-glass transitions widely observed in pyrochlore- 
based frustrated materials, we investigate effects of coupling to local lattice distortions in a bond- 
disordered antiferromagnet on the pyrochlore lattice by extensive Monte Carlo simulations. We 
show that the spin-glass transition temperature Tf is largely enhanced by the spin-lattice coupling, 
and furthermore, becomes almost independent of A in a wide range of the disorder strength A. The 
critical property of the spin glass transition is indistinguishable from that of the canonical Heisenberg 
spin glass in the entire range of A. These peculiar behaviors are ascribed to a modification of the 
degenerate manifold from continuous to semidiscrete one by the spin-lattice coupling. 

PACS numbers: 75.10.Hk, 75.50.Lk, 75.10.Nr 



In the last few decades, increasing attention has been 
devoted to low-temperature (T) behavior of geometrically 
frustrated magnets [l|. Frustration suppresses conven- 
tional long-range ordering such as Neel ordering down to 
much lower T compared to the interaction energy scale, 
opening the possibility of alternative low-T phases. The 
spin glass (SG), in which spins are frozen randomly, is 
one of such possibilities widely observed in geometrically 
frustrated materials However, it is unclear so far 

how the nature of SG is different from the canonical one 
driven solely by randomness 

Antiferromagnets on a pyrochlore lattice (the inset 
of Fig. [1]) are typical examples of geometrically frus- 
trated spin systems. When considering classical Heisen- 
berg spins with nearest-neighbor exchange interactions, 
no long-range ordering occurs down to zero T, and 
the ground state has continuous macroscopic degener- 
acy m E3- Recently, the effect of randomness in the 
exchange interactions was studied upon this degenerate 
manifold ]j]-14|. It was shown that the randomness im- 
mediately lifts the degeneracy, inducing a SG transition: 
The transition temperature Tf is proportional to the dis- 
order strength A. This gives a clue to explain why SG is 
prevailing in geometrically frustrated materials. 

However, several characteristics of the SG still remain 
puzzling. One of the surprising aspects is that, in many 
pyrochlore-based magnets, Tf appears to be almost in- 
dependent of the strength of disorder A. For example, 
for a typical frustrated SG compound Y2M02O7 with 
T f ~ 22K [HHIl, random substitution of Y 2 + by La 2+ 
does not change Tf for the La concentration up to 50%, 
despite a substantial change of the Curie- Weiss temper- 
ature Oqw Hoi- Similar behavior was observed also in 
spinel oxides (Zni^Cd^C^C^ for x < 0.1, in which SG 
emerges after the antiferromagnetic spin-lattice order at 
x = is destroyed by small Cd substitution [21]. An- 
other distinctive aspect is that Tf is much higher than 
that theoretically expected for a moderate strength of 
disorder A; e.g., a numerical estimate of Tf/J ~ 0.01 for 
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FIG. 1: (color online). A-T phase diagrams obtained by MC 
simulation with (a) J™ op = 0.0 and (b) 0.075 for b = 0.2. 
The nematic (T c ), antiferromagnetic (Tn), and SG transition 
temperatures (Tf) are denoted by squares, triangles, and cir- 
cles, respectively. In (a), Tf coincides with T c for A > 0.3, 
suggesting a multicritical point at A ~ 0.3 (~ b). The cross in 
(a) denotes Tf for b — and A = 0.1 14]. The inset shows a 
16-site cubic unit cell of the pyrochlore lattice. The open and 
filled circles in the inset denote two nonequivalent sites with 
opposite spins in the Neel order. The Neel-SG phase bound- 
ary in (b) was determined by the inflection point of m s 2 (A) 
curve for L — 4 [see the inset of Fig. 2(f)]. 



A/ J = 0.1 is about 20-30 times smaller than the exper- 

rs indicate that Tf 



imental value 

is not set by A, and suggest that some important factor 
is missing in the previous theories [nun. 

A candidate is the magnetoelastic coupling. Impor- 
tance of local lattice distortions has been pointed out for 
Y2M02O7 by various microscopic probes [22l426j. They 
are crucial also in (Zni-^Cd^C^O^ as obviously seen 
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in the spin-lattice ordering at x = [2l| . Theoretically, it 
was shown that the bond randomness destroys the spin- 



lattice ordering and induces a SG state [27j. However, 
the argument was limited to a uniform global lattice dis- 
tortion, and Tf is deduced to behave similarly to the case 
in the absence of the spin-lattice coupling, i.e., Tf ex A. 
Hence, the puzzling behaviors of Tf still remain unre- 
solved. 

In this Letter, we investigate the effect of coupling to 
local lattice distortions on the SG transition in a bond- 
disordered classical Heisenberg antiferromagnet on the 
pyrochlore lattice. We show that Tf is largely enhanced 
by the spin-lattice coupling 6, and pinned at T c± 6, i.e., 
almost independent of A in a wide regime of A. We clar- 
ify the mechanism of the peculiar behaviors by consider- 
ing how the spin-lattice coupling affects the degenerate 
manifold in the model. 



We take as a starting point the Hamiltonian 



-> -> K 
Jij (1 - apij) Si • Sj + —p% 



(1) 



where the sum runs over the nearest neighbor bonds of 
the pyrochlore lattice, and pij is the change in distance 
between neighboring Heisenberg spins Si and Sj, rela- 
tive to the equilibrium lattice constant. The model de- 
scribes both the exchange randomness in Jij induced by 
static bond disorder, and magnetoelastic coupling to lo- 
cal lattice distortions . Here we introduce a uniformly- 
distributed randomness as Jij G [J — A, J + A] with 
< A < J. The last term in Eq. (pQ) represents the elas- 
tic energy of lattice distortions. Hereafter, all the energy 
scales including T are measured in units of J. 

In general, the lattice distortions pij depend on each 
other; that is, they have a cooperative aspect, which may 
lead to structural transition concomitant with some mag- 
netic ordering. When the cooperative aspect is ignored, 
the model (pQ) is much simplified by integrating out pij. 
The resultant spin-only Hamiltonian is given by 
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Jij Si • Sj bij (^Si • Sj 



(2) 



where is the biquadratic coupling given by 

J 2 J a 2 /2K > 0, also being bond-disordered variables. 
Hereafter we use b = a 2 /2K to measure the strength of 
the spin-lattice coupling. We consider the model (|2]) to 
unveil intrinsic effects of b. The result will be discussed 
for Mo pyrochlores that show no spin-lattice ordering. 
To discuss the case in which a spin-lattice ordering ap- 
pears as in (Zni-^Cd^C^O^ we take the cooperative 
aspect into account in a form of effective spin-spin inter- 
actions, following the results in previous theoretical stud- 
ies [111. Among many possible contributions, we focus on 
the simplest one from anti-correlations between neighbor- 
ing bond distortions by introducing a next-nearest neigh- 
bor term H 

T-coop 
J 2 



coop 

= ^2 COOP E(( M -)) Si • Sj to Eq. © with 
> (the sum is over the second neighbor pairs; see 



the inset of Fig. [T]) (39|. A similar term is derived from 
the coupling to the site phonon [29|. 

When b and J2 00p are both zero, which we call the 
bilinear limit, the model reduces to the one previously 
studied Ill4l4j. For a finite spin-lattice coupling b > 
but J2° op = 0, the present model at A = exhibits 
a weak first-order nematic transition at T c ~ 6, below 
which spins select a common axis without selecting their 
directions on it. The ground state is identified by a set 
of spin-ice type local constraints that two of four spins 
are opposite to the other two on every tetrahedron. This 
leaves the system magnetically disordered down to T = 
with the semi- discrete degeneracy [30|. J^ 00 ^ lifts the 
degeneracy and stabilizes a spin-lattice (Neel) order over 
the nematic phase. The antiferromagnetic J^ 00 ^ induces 
the q = four-sublattice Neel order (3l[, schematically 
shown in the inset of Fig. HJ Our interest here is how the 
SG transition is induced by A in the competition with 
the nematic and Neel orderings. 

To address the issue, we employ an extension of 
the loop algorithm which enables an ergodic sampling 
over the semidiscrete degenerate manifold at low T <C 
b [HI. We also adopt the exchange MC method [HI 
and the overrelaxation update [34|. We consider pe- 
riodic systems of cubic geometry with L 3 unit cells 
with totally N s = 16T 3 spins. To identify the SG, 
nematic, and antiferromagnetic transitions, we calcu- 
late the SG susceptibility Xsg, nematic order param- 
eter Q 2 , sublattice magnetization ra s , and specific 
heat C. xsg is given by N s q^ A , where #| A is the 
Edwards-Anderson order parameter defined by g| A = 
^((E^yJE^SrX^A [3l. Here (• • -) T 
denotes a thermal average and (• ■ -)a a random aver- 
age over the interaction sets; the upper suffixes a and j3 
denote two independent replicas of the system with the 
same interaction set. The nematic order parameter Q 2 , 
which measures the spin collinearity, is defined as Q 2 = 
2^s" 2 ((E^{(^ ' Sj) 2 - 1/3})t)a- The sublattice mag- 
netization is defined as ra s = \(J2iei Si)\ 2 /N 2 ) 1 ^ 2 , 
where I labels the four sublattices of the pyrochlore lat- 
tice [3lj . The specific heat C is calculated by C = 

N~ 1 T~ 1 (('H 2 }t - (W)t)a. An data shown below are 
averaged over a number of interaction sets varying from 
2000 to 100, and typical MC steps for measurement vary 
from 10 4 to 10 T depending on L and A. 

We start from the result in the absence of J2 00p . A typ- 
ical phase diagram is presented in Fig.[]Ja) with b = 0.2. 
The system exhibits a nematic transition at T c ^ 6, 
but remains magnetically disordered down to zero T at 
A = 0. By introducing the disorder A, the SG transition 
appears at a finite Tf. In the weakly-disordered region, 
i.e., A < 6, Tf is roughly proportional to A; we call this 
regime the linear regime. A remarkable point is that Tf is 
largely enhanced compared to that in the bilinear limit: 
the enhancement factor is, e.g., about 5-10. At A c± 6, 
Tf appears to merge into T c with showing multicritical 
behavior. For larger A, Tf (= T c ) becomes nearly inde- 
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FIG. 2: (color online). Temperature dependences of the ne- 
matic order parameter Q 2 , specific heat C, SG susceptibility 
Xsg, and sublattice magnetization m s for J^ 00 * 3 = [(a)- (d)] 
and for J 2 coop = 0.075 [(e) and (f)] at b = 0.2. The insets in 
(b) and (d) show the scaling collapses of xsg- The inset in 
(f) presents ra 2 (A) at T = 0.25, 0.2, and 0.1 for L = 4. 



pendent of A, which we call the plateau regime. This is 
in sharp contrast to the previously-reported SG behav- 
ior, Tf ex A [l2l-[l4). These peculiar dependences of Tf on 
A are universally obtained for a wide range of < b < J. 

Figures [2(a) and [2(b) show typical MC data in the 
linear regime A < b. As seen in Fig. [2(a), Q 2 exhibits a 
steep rise at T c± 0.22, signaling the nematic transition. 
At the same T, the specific heat C exhibits a sharp peak. 
By extrapolating the peak position to the bulk limit, we 
obtain T c = 0.2245(15) at A = 0.1. On the other hand, 
Xsg appears to diverge at a lower T, indicating the SG 
transition. We obtained T f = 0.102(14), 7 = 2.24(75), 
and v = 1.16(18) by the finite-size scaling analysis as 
demonstrated in the inset of Fig. [2(b): Here, we assume 
Xsg = L^< v f{L x l v t), where t = (T — T f )/T f , v and 7 are 
the critical exponents for the correlation length and xsg, 
respectively. The critical exponents are consistent with 
those in the bilinear limit 6 = 12|-14] as well as of the 
canonical SG [§[ within the error bars, but surprisingly, 
Tf is strongly enhanced by a factor of 5-10 from Tf = 0.01- 
0.02 at 6 = [TMl- 



What happens in this linear regime is the following. 
The spin collinearity emergent below T c enforces spins 
to satisfy the spin-ice type local constraints and form 
locally-correlated collinear objects. There, the system 
bears a semidiscrete degenerate manifold with a multival- 
ley energy landscape, as in the case of A = 0. Therefore, 
thermal fluctuations are strongly constrained compared 
to the bilinear case in which the ground-state manifold is 
continuously-connected. At the same time, the spin-spin 
correlations are much enhanced to exhibit quasi-long- 
range behavior below T c [3(j • These are the origin of the 
striking enhancement of Tf demonstrated in Fig. QJa). 

With further increasing A, Tf reaches T c at A ~ b but 
never exceeds T c [40]. This is understood by the above 
physical picture. Typical data in the plateau regime are 



shown in Figs. [2(c) and [2(d). Following the above way, 
we obtained T c = 0.256(3) and T f = 0.252(4) at A = 0.5. 
The scaling analysis of xsg is compatible with the second 
order transition with critical exponents 7 = 1.2(5) and 
v = 0.67(16) [see the inset of Fig. [2(d)], which are also 
consistent with those in the bilinear limit [4l|. These 
results indicate that the SG transition coincides with the 
onset of spin collinearity, and the concomitant transition 
looks second order in this plateau regime. The critical 
behavior is indistinguishable from the conventional SG [H, 
12- TiJ: It is noteworthy that the peak of C is suppressed 
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and broadened [Fig. [2(c)], as seen in the canonical SG 
transition jsj. 

Next, we examine the effect of the cooperative coupling 
. A typical phase diagram is shown in Fig. H(b) with 
= 0.075 and b = 0.2. At this value of J 2 coop , the 
nematic phase in the linear regime is taken over by the 
q = Neel-ordered phase. Figures [2(e) and [2(f) indicate 
a first-order transition at Tn = 0.291(3) for A = 0.2. On 
the contrary, the SG in the plateau regime remains ro- 
bust; we confirmed that the MC results in this regime are 
essentially the same as at J^ 00 ^ = shown in Figs. [2(c) 
and [2(d). Thus, J^ 0013 lifts the residual semi-discrete de- 
generacy in the nematic phase, hardly affecting the SG in 
the plateau regime. This results in a bicritical-like phase 
competition between the Neel and SG phases. 

The robust plateau behavior of Tf at a largely enhanced 
value gives an explanation for the puzzling behaviors in 
the pyrochlore-based antiferromagnets mentioned in the 
introduction. For example, in (La^Yi-a^Mc^Or, many 
experiments suggest a substantial bon d disor der even in 

Suppose 
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the stoichiometric case at x = 
that the system inevitably includes a substantial disorder 
even in the best-quality sample ever made, as suggested 
by the microscopic probes 22, 23|, and is already in the 
plateau regime, the system does not pass through the ne- 
matic phase and Tf can be large and remain almost con- 
stant against the additional disorder by La substitution 
x. In the case of (Zni-^Cd^C^C^, the results at a fi- 
nite J 2 give a qualitatively reasonable explanation for 
the phase competition between the spin-lattice ordered 
phase and SG phase as well as the peculiar behavior of 
Tf. For further qualitative comparison, more sophisti- 
cated modeling is indispensable for farther neighbor ex- 
changes, substitution effect, and spin-lattice couplings. 

Finally, we make a further comparison with 
experiments by the uniform magnetic susceptibil- 
ity. The susceptibility, calculated by xo = 

is shown for various values 
in Fig. [3(a). The results indicate 
that xo show a Curie- Weiss-like behavior at much higher 
T than Tf with a A-dependent Curie- Weiss temperature 
#cw- As shown in Fig. [3(b), |#cw| decreases significantly 
with increasing A, whereas Tf stays almost constant in 
the plateau regime. This constrasting behavior of Oqw 
and Tf is robust for finite J^ 00 ^. The result qualita- 
tively agrees with the peculiar SG behavior observed in 
(La^Yi-^JVk^Oy [20[. Note that similar behaviors of 



of A with J 2 coop 




FIG. 3: (color online), (a) Inverse magnetic susceptibility 
X^ 1 at J2° op = for L — 3. The lines denote Curie- Weiss 
fittings for 0.6 < T < 0.9. (b) A dependences of 6>cw and T f . 
The circles and squares denote the results for J^ 0015 = and 
0.075, respectively. 



Tf and Oqw were reported in i?2Mc>2 07 for a change of 
R cations, such as Gd, Tb, and Dy 3(1 137|. 

To summarize, through the numerical analysis of the 
bond-disordered pyrochlore antiferromagnet, we have re- 
vealed that the SG transition temperature Tf is largely 
enhanced by the spin-lattice coupling 6, and becomes sat- 
urated at a temperature set by b. The results reproduce 
the puzzling characteristics of the spin-glass transition 
observed in the pyrochlore-based geometrically frustrated 
materials. Unfortunately, the value of b is not available 
for the compounds that we are concerned here [42|, but 
we believe that our results stimulate experiments to esti- 
mate of 6, e.g., from high-T nonlinear magnetic suscepti- 
bility (38]. Since the peculiar behaviors originate in the 
simple local physics, we believe that they are common to 
systems coupled to local lattice distortions. 
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